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Abstract 

We study theoretically the quantum dynamics of an electron in the singly- 
ionized double-donor structure in the semiconductor host under the influence 
of laser pulses whose frequencies are close to structure resonant frequencies. 
This system can be used as a charge qubit where the logical states are defined 
by the lowest two energy states of the remaining valence electron localized 
around one or another donor. The quantum operations are performed via 
resonant or Raman-like transitions between the localized (qubit) states and 
the excited states delocalized over the structure, combined with phase shifts 
between qubit states generated by voltage pulses. The possibility of realization 
of arbitrary single-qubit rotations is demonstrated. 
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I. INTRODUCTION 

In view of recent progress in the development of the controlled-donor implanta- 
tion techniques Kane's paradigm of the solid-state quantum computations [2] 
has gained new insights. The alternative schemes using the orbital [3|, H] or spin 
[5l E] degrees of freedom of the donor-bounded electrons to encode the quantum in- 
formation instead of the nuclear donor spins, have been proposed. Besides, several 
refinements of the original proposal concerning the initialization [7J and read-out 
[El [9l [TOl [11] as well as the information transfer through the quantum networks 
[ni [131 [H] have been made. 

In particular, the pair of donors sharing an electron has been considered as very 
promising candidate for the solid-state qubit embodiment [5]. The qubit is pre- 
sented by the electron orbital states positioned at the different donors. There are 
two main driving mechanisms for the coherent electron evolution defining a quan- 
tum operation on such a qubit. First uses the electric fields through an application 
of the adiabatically switched voltages to the surface gates placed above the donor 
structure to modify the confinement potential, thus varying the electron tunneling 
rates between neighboring donors [3l [151 [IS] • The desired state of the qubit is real- 
ized by an appropriate choice of the voltage parameters. The second scheme relies 
upon the optical dipole transitions between the size-quantized one-electron levels 
induced by the resonant pulses [TTl [TH|. This technique requires one to irradiate 
the donors encapsulated into near-surface semiconductor layers by laser pulses with 
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frequencies lying in terahertz range. Such kind of radiation sources is currently 
under extensive exploration. To all appearances, the quantum cascade lasers [IS] 
in which photon generation takes place via electronic intersubband transitions in 
semiconductor heterostructures, will be able in near future to cover this frequency 
range. As it was shown, an arbitrary single-qubit operation can be achieved with 
two simultaneously switched pulses connecting the qubit states via the excitation 
of the state delocalized over the structure. The latter scheme is likely to be more 
preferable than the former due to its higher selectivity and lower field intensity. 
The successful implementation of the quantum operations, however, requires one to 
provide the high precision in the durations, the frequencies, the polarizations and 
the strengths of the pulses together with reliable control over the delay time and 
the phase difference between the pulses. Besides, the use of the intermediate state 
lying in the neighborhood of the continuum introduces the decoherence caused by 
the ionization and the spontaneous emission from this state. 

Here we propose the way to overcome the difficulties inherent to the electron 
charge manipulations by optical means. We show that the off-resonant laser pulses 
combined with electrostatic fields, can be used to generate an arbitrary rotation of 
the qubit-state vector as well as to drive an electron between the remote donors. 
It is essential that this method is based on the Raman-like transitions between the 
localized electron states of an effective molecular ion where the excited (delocalized 
over the double-donor system) states irrespective of their number are used as the 
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transport channels. As we shall see, the coherent electron dynamics is described 
by the simple analytical model. In order to check analytical results, we have per- 
formed numerical simulations on single-electron dynamics in the effective hydrogen 
molecular ion. This system is very popular object for modelling electronic proper- 
ties of real donor structures (see, e.g., Refs. [18] and [20]). The eigenenergies and 
eigenfunctions of molecular ion were found from the stationary Schrodinger equation 
and then used to calculate the set of eigenfrequencies and the matrix elements of 
optical dipole transition between eigenstates of hydrogen molecular ion that define 
the field-structure interaction strengths. The time- dependent Schrodinger equation 
generating the coherent electron evolution in the external fields was integrated nu- 
merically for different values of structure and field parameters. By comparison of 
analytical expressions of the Rabi frequencies for both resonant and off-resonant 
cases, we found them to be in good agreement with those extracted from numeri- 
cal solutions. This result confirms the validity of approximations made during the 
analytical treatment. The pulse and structure parameters needed for those types of 
quantum evolution may be evaluated from the results obtained in this study. 

The paper is organized as follows. In Section II we describe general model of 
coherent electron dynamics in double-donor structure under the influence of near- 
resonant electromagnetic pulse. The Raman-like off-resonant electronic transitions 
are studied analytically in Section III. Besides, in this Section the possible realiza- 
tion of several basic single-qubit quantum operations is proposed. The results of the 
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numerical study of electron dynamics in hydrogen molecular ion are given in Sec- 
tion IV. The advantages of the off-resonant driving scheme in comparison with the 
resonant scheme proposed earlier as well as some questions concerned with further 
exploration of optically-driven donor-based charge qubits, are discussed through the 
Section V. We conclude our paper by Section VI. 

II. MODEL 

We begin with the description of the one-electron double-donor (DD) structure 
(Fig. 1 (a)). Let the donors A and B be placed on the axis z from each other 
at the distance R (hereafter - internuclear distance) large enough to consider their 
ground orbital states to be isolated. Due to this fact those states may be used as 
the qubit states |0) and |1) (for example, if an electron is localized on the donor 
A{B), the qubit is in the state |0) (|1))). The energy difference A — £i — may 
be introduced due to the structural asymmetry caused by the fabrication process 
and/or via the bias voltages generated by the surface gate Va- The coupling be- 
tween the excited states of individual donors through the electron tunneling gives 
rise to the forming of hybridized states delocalized over the DD structure. However, 
if it! » (a^ is the effective Bohr radius of host material) the low-lying excited 
states are hybridized weakly and do not participate the two-donor dynamics. We 
shall be interested therefore only in the excited states of individual donors whose 
wave functions considerably overlap. The single-electron spectrum of the DD struc- 
ture is presented by the sequence of the states {\k)}^^^ which for A = are the 



5 



doublets composed of the symmetric and antisymmetric superpositions of isolated 
donor states. If A 7^ 0, the spectrum is assumed to be much more complex, as 
one can see from Fig. 1 (b). Taking into account the hydrogen-hke spectrum of the 
isolated donors we expect the excited states close to the edge of the potential barrier 
separating the donors to have the quasi-continuous energy distribution. 

In what follows we shall study the one-electron quantum dynamics involving the 
localized (qubit) states, |0) and and the states {\k)}^^^ delocalized over the 
structure. Our aim is to choose the field and structure parameters so that to drive 
an initial qubit state |^ {to)) = ao |0)-|-/?o |1) = {do, Po)^ into the final state |^ {t)) = 
a\0) + P\l) — {a, (5)^ with the desired coefficients a and (3. 

In the absence of an external field the DD structure is characterized by the sta- 
tionary Hamiltonian Hq with the eigenstates {\k)} and the eigenenergies {sk}'- 

Ho \k) = Sk \k) . (1) 
The eigenstates {\k)} form the complete orthonormal set so that 

j:\k){k\^i. (2) 

k 

In the presence of the electromagnetic field the system Hamiltonian reads 

H = Ho-eE {t) r, (3) 

where e is the electron charge, E {t) is the field strength, r is the radius- vector of an 
electron. With the help of Eqs. (1) and (2) we rewrite the Eq. (3) in terms of the 
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projection operators: 

^ = (E I^) (^i) H \m) (ml) = ^ £fe (^| + E (t) ^ d,^ |^) (m| , (4) 

\ fe J \m J k k,m 

where dkm — {k\ — er \m) is the matrix element of optical dipole transition between 
the states and \m). The state vector of the system may be presented in the form 

l*W) = EcnWe-"*|n) (5) 

n 

and is governed by the non-stationary Schrodinger equation 

.M = /f|*W>. (6) 

with the initial condition |^ (to)) = olq |0) + /3o |1) (hereafter h = 1). 

Inserting Eqs. (4) and (5) into Eq. (6), we arrive at the set of linear differential 

equations for the probability amplitudes c„ (t) . We shall only examine the transitions 

between the states |0) and |1) and the states {\k)}^^^: 
ico = E(i)EdofeCfee-^'^o'=* 

k 

ici = E(i)EdifcCfee-^'^i'=* (7) 

ick = E (t) (dSfcCoe-o'=* + d^,cie--*) , k G {A;},,, , 

where a;o(i)fc = £k - £o(i)- 

Let the electromagnetic field imposed on the structure to have (in the dipole 
approximation) the form of two phase-locked pulses 

E {t) = Eo {t) cos {ujot + cpo) + El {t) cos {uit + cpi) , (8) 

where the pulse envelopes Eq (t) — Eq/o (t) , Ei (t) — Ei/i (t) are the slowly- varying 
(compared to optical frequencies) time-dependent functions, a;o,i are the pulse fre- 



quencies, and (/7o,i are the pulse phases. We require both pulses to be in the two- 
photon resonance with the DD structure, i.e., Eq + luq = Si + lui or, alternatively. 



^Qk — ^Ik 



(9) 



where 5o(i)fe = a;o(i) — (^o(i)fe is the detuning of the pulse frequency a;o(i) from the 
resonant frequency a;o(i)jt. Note, that we use two pulses only if A 7^ 0. In the 
symmetric structure, where A = 0, the single pulse E (t) = Eg (t) cos {ujQt + ipo) is 
enough. 

Making use of the rotating-wave approximation (that implies A -C ci;o(i) and 
\5k\ <S <^o(i)), we obtain from Eq. (7) the following set: 



k 
k 

iCk^\\lk{t)+M)(^''^' 



r,iSl,t 



^^lk it) e-'^' + Kk it)] cie-^^'=*, k e , 



(10) 



where Ao(i)fc (t) = Ao(i)fc/o(i) (t) e"^°w , /Xo(i)fc {t) = i^o{i)kfo{i) (t) e^'^od) , Ao(i)fc = do(i)fcEo(i)/2, 
Aio(i)fc = di(o)fcEo(i)/2 and the identities a;o(i) - u;^o)k = (^i(o)fc =F A are used. 

Eq. (10) describes the dynamical process involving many three- level excitation 
schemes that act in parallel. Each of them is characterized by the set of parameters 
A, Ao(i)fe, A*o(i)fe, and 5k, where k e {k}g^^- We shall suppose the values of Ao(i)fc and 
A*o(i)fe to be of the same order. Depending on the ratios between these parameters, 
k-th excitation scheme may be classified in the following way. First we consider the 
case of small detunings. If the coupling coefficients of the optical dipole transitions 
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Ao(i)fc and the detunings 5k satisfy the inequahty 

\h\ , A < |Aofc| , |Aifc| , (11) 

the three-level scheme works in the resonant symmetric regime. Instead, the appli- 
cability of the resonant asymmetric scheme [TTj is provided by the condition 

141 < |Aofc|,|Aifc| < A. (12) 

We see that the asymmetry/ symmetry of the structure isn't defined by the pres- 
ence/absence of the energy difference A only but by the ratio between A and the 
coupling coefficients |Aoyfc| , |Aijt| as well. In other words, a driven DD structure can 
be treated (relative to the fc-th transition scheme) as symmetric if the infiuence of 
the parameter A, introducing a "static" asymmetry, is compensated by an appro- 
priate value of the field strength defined from (11). In this case only one external 
pulse is sufficient to excite both transitions [2T] . 

Next we shall examine the opposite case where the states |0) and |1) are connected 
through the off-resonant transitions involving a collection of the excited states lying 
at the edge of the potential barrier. Two situations are possible again: 

A < |Aofe|,|Aife| < (13) 

and 

|Aofc| , lAifcl < 141 , A. (14) 

The first of these inequalities corresponds to the off-resonant symmetric excitation 
scheme. This situation was studied in Refs. [221 123] for the double-dot structures 
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with three and four levels. As it was shown, the set of single-qubit operations 
produced by such one-electron dynamics is incomplete, since in order to realize an 
arbitrary rotation of the qubit-state vector, the structure symmetry must be broken. 
In what follows, our attention will be focused on both the symmetric case and the 
off-resonant asymmetric case for which the conditions (14) are satisfied and each 
pulse drives the transitions between only one of localized state and the transport 
states {|fc)}e2;c- This implies also that the values of A and 6k must be rather different 
from each other for all k to prevent the single-donor resonant dynamics. 

III. THE OFF-RESONANT DYNAMICS 
In this Section we consider two cases of the off-resonant electronic quantum dy- 
namics in DD structure, i.e., the off-resonant symmetric dynamics and the off- 
resonant asymmetric dynamics. We begin with the second one, initially specifying 
the conditions imposed on pulse and structure parameters: 



0(l)fc 



/^0(l)fc 



< A < 
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S/iO(l)/c 



A 



0(1)A: 



/^0(l)fe 



(15) 



(The case of large A is out of scope of this study; see Sec. IV D for the reasons). 

Strictly speaking, the rotating-wave approximation used in derivation of Eq. (10) 
and based on the averaging of Eq. (7) over time interval [r — tt/uq, t + it/uq] (RWA 
1) [211, cannot be applied to high-lying states for which the detunings 6k become 
comparable to the pulse frequencies. Let us define effective "maximal" detuning 
6max < ^0 that corresponds to the upper bound of energy range, where RWA 1 is 
still valid, and average Eq. (10) over time interval [r — 7r/(5max, t + vr/(5inax] (RWA 
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2). This step is justified by first of tlie inequalities (15). After some manipulations, 
we obtain two equations for time-averaged functions Cq and Ci: 

^co = E ^ {|Aofc + /iifc exp {-iA ■ t)f cq + [Xok + /^ife exp {-iA ■ t)] [Ai^ + fiok exp {iA ■ t)]* ci}, 
ici = J2j^ ||Aifc + /iofcexp {iA ■ t)f ci + [Aofc + /ii/t exp {-iA ■ t)]* [Ai^ + /iofcexp {iA ■ t)] cq}. 

(16) 

Again, careful consideration of averaging procedure requires us to handle only with 
the detunings that are sufficiently small in comparison with (5max- It is easy to see 
that the equations (16) do not contain the fast-oscillating terms Cfc^o,i- Therefore, 
these equations describe effective two-state dynamics involving only the localized 
states |0) and |1). From third equation of the set (10) we conclude that the functions 
Cfc^o,i have the amplitudes ~ |AoA:/5fc| << 1 and oscillate at the frequencies ~ \Sk\. 
The reduction of set like Eq. (10) to the two equations for slowly-varying probability 
amplitudes Cq and Ci is known in atomic optics as the adiabatic elimination procedure 
[23] - To proceed further, let us average the equations (16) over the time interval 
[r — tt/A, r + tt/A] (RWA 3) on which all of the time- dependent functions except 
exp {±iAt) may be replaced by their mean values so that after integration we arrive 
at the following set of equations: 

ico = T.j-{ [|Aofc {t)f + l/iifc (t)l^l Co + Aofc (t) X*ik (t) cij, 



(17) 

^ci = E i {Kk (t) Aifc {t) Co + [lAifc (t)l' + l/iofc (t)l'] ci}. 



k 

E 

k 

Now we return to the symmetric driving scheme that is simpler than the asym- 



metric one. By setting A = and Ei = 0, we transform the set of Eq. (10) 
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into 

ico = E ^ofe (t) hk 

k 

ici^j:Xik{t)bk (18) 

k 

where bk — exp (iSkt) , k e {^}ea;c ^® redefine the couphng coefficients as 



Ao(i)fc (t) = Ao(i)fe/o (t) e"^°, Ao(i)fc = do(i)feEo/2 . The inequahties (13) allow one to 
apply the adiabatic elimination procedure to the intermediate levels {\k)}^^^: 

bk ^0,bk^ [XI, {t) Co + A*, (t) c,]/5k, k e {fc},,, . (19) 

The equations for two remaining probability amplitudes Cq and Ci in the matrix 
form for both asymmetric and symmetric dynamics read 

co)_( Ao{t) A,{t)\(co\ . . 

where for symmetric case one has Aq (t) = Aq/q (i), Ai (i) = Ai/Q(i), A2 (t) = 
A2/2 (t) and 



^o = J2 |AofclY ^fe, Ai = J2 |AifelY 4, A2 = J2^okX*ik/^k: (21) 

k k k 

whereas for asymmetric case the functions entering into Eq. (20) are given by expres- 



k 



sions Ao(t) = E \Xok{t)\ +\nik{t)\ /4, Ai (t) = E |Aife(t)| +|/xofc(^)| /4, 



k 



|2 , I /,m2 



and A2 (t) — E Aoifc (t) XI, {t)/5k. It allows us to treat both cases in the same man- 

k 

ner, however, keeping in mind that for symmetric structure Aq (t) = Ai (t), since 
Aofc = ±Aife for the state |A;) whose electronic wave function possesses even/odd par- 
ity relative to the center of DD structure, that does not take place for asymmetric 
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structure. Below we present the general solution of Eq. (20) that will describe both 
types of off-resonant electron dynamics. 

The eigenstates and the eigenenergies of the matrix in right-hand side of Eq. (20) 
may be written as 



1+) = e^^'-st^^Wl cos [0 (t)/2] |0) + sin [0 (t)/2] |1) 
I-) = e'^'-s[A2W] sin [0 {t)/2] |0) - cos [0 \t)/2] |1) 



(22) 



and 

E± {t) = [Ao {t) + Ai {t)]/2 ±n{t), (23) 

respectively, where 

cos [0 it)] = [Ao (t) - A^ {t)]/2n (t), sin [0 (t)] = jA^ {t)\/Q (t), (24) 

and 



n{t)^^ [Ao (t) - Ai {t)f/i + \A2{t)f (25) 
is the instantaneous Rabi frequency. Using the unitary transformation 

n(A- ( e'^-'st^^W] cos [0 (t)/2] e^'^-'st^^W] gin [0 (t)/2] \ , . 

^^^>-[ sin [0 (t) /2] - cos [0 (t) /2] ) 

we represent the state vector in the instantaneous basis {|-|-) , |— )} as 

1$ (t)) = a+ (t) 1+) + a_ (t) I-) , 1^ (t)) = D (t) 1$ (t)) (27) 
and rewrite Eq. (20) in the new basis as 
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where 



e{t) 



[Ao (t) - Ai (t) 


IA2WI- 




[Ao (i) - A^ (t)] 


[Ao (t) 


-Ai it)]'/A+\A,it)\' 



(29) 



Here we restrict our interest by the diagonal evolution followed from the choice of 
the system parameters for which © (t) <^ E± (t) . In this case the solution of Eq. 
(28) is straightforward: 

a+ {t) = a+ (to) exp 
a_ (t) = a_ (to) exp 

With the help of equations (26), (27), and (30) we may write down the expression 
for the evolution matrix of the qubit-state vector |^ (t)) in the laboratory frame: 
\^{t))^U\^{to)), 



-i J E+ {t') dt' 

to 

-i J E_ {t') dt' 

to 



(30) 



U 



-isot 







-ieit 



Dif) 



^-iieot+fKit)] 



«00 



-i J E+{t')dt' 

e *o 



\ 





t 

-i J E-(t')dt' 

e *o 



D"^ (to) 



(31) 



Uioe Uiie j 



where 



Uqq = = e-^^(*) cos [e {t) /2] cos [O {to) /2] + e^^(*) sin [0 {t) /2] sin [O (to) /2] , 



= e^arg[A2(t)] |g-in(t) [Q (^)/2] sin [O (io)/2] - e^^W sin [0 (i)/2] cos [O (to)/2]} , 

(32) 

and 

t t 
Q (i) = 1 Q (i') dt', (t) = y [Ao (f) + Ai (i')]/2di'. (33) 

to to 

The expressions (31) - (33) describe the effective two-level dynamics that corresponds 
to the continuous evolution of the qubit state vector on the Bloch sphere. In the 
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remainder of this section we show how to choose the pulse and structure parameters 
in order to reahze the most important single-qubit gates. We illustrate the qubit 
state engineering by considering a particular case of the driving pulses sharing the 
same time dependence, i.e., /o {t) — fi {t) = f {t). The condition 6 = is then 
satisfied and the components of the evolution matrix (31) take the form 



1*00 = 1*11 — cos 0,{t) — i cos (©o) sin fl (t) 
= -ulo = -ie^^g(^2) sin (Oo) sin iCl (t) 



(34) 



where 6o = arcsin 



A2|/J(Ao-Ai)74+|A2 



. The dynamics described by the 



equations (34) is sufficient to generate an arbitrary single-qubit rotation on the 
Bloch sphere. For example, the quantum operations such as NOT {ax): {ckq, (5q) — > 
(/3o, cto) ; PHASE (cr^): (ao,/3o) ^ (ao, -/3o) ; and Hadamard (if): {ao,Po) 
(cKo + Po) /V^, (cKo ~ Po) /V^ can be realized (up to the common phase) given 
the following choices of the pulse - structure parameters: 



n (T) = 7r/2 + TTk, TA = 2t:1, arg (A2) = 27rm 



(35) 



and 



©0 ((^x) = 7r/2 + 27rn, Gq (ctz) = vrn, 0o (H) = n/A + Tin, 



(36) 



respectively. Here k, 1,171, and n are the integers and T is the pulse duration. Of 
course, this is not a unique parameter choice to attain the above quantum operations. 
Careful estimation of the operation times requires the detailed knowledge of the 
energy spectrum of DD structure and the values of Ao(i)fe. 

IV. A HYDROGEN MOLECULAR ION: NUMERICAL STUDY 
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We visualize our results with the help of model, where the energies and wave 
functions of an electron, bounded in the DD structure, are approximated by the 
eigenenergies and eigenfunctions of effective hydrogen molecular ion H^, for which 
the effective Bohr radius = e(mo/m*)aB and the effective Rydberg energy Ry* = 
[m* /{e^mQ)]Ry contain the information about real solid-state environment via the 
electronic effective mass m* and the dielectric constant e [Ry = 13.6057 eV is 
Rydberg energy, as = 0.5292 x 10^^° m is the Bohr radius, and mo = 9.1094 x 10^'^^ 
kg is the free electron mass). In practice, the single- valley approximation was applied 
to study the phonon-induced decoherence of electron in Si:P^ structure [20] and 
the effect of surface gate on single phosphorous donor in the silicon [28], as well 
as to calculate the Rabi frequency of resonant electron transfer between localized 
states in SiiPj' structure in three-level approximation [TH]. The comparison of the 
results of calculations, where intervalley interference is taken into account, with those 
neglecting such effects, demonstrates qualitative agreement between them (especially 
in presence of an external field, see [IS]). In our study, we are mostly interested 
in validation of the approximations, made during the paper to reveal the specific 
features of resonant and off-resonant electron dynamics in a multilevel structure. 
We expect, that the dynamical analysis developed above and confirmed numerically 
below for effective single- valley-approximated structure, can be applied to real solid- 
state DD structures as well. 

In what follows we show how to find the eigenenergies and eigenfunctions of 
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hydrogen molecular ion (HMI) without LCAO approximation. We calculate the 
matrix elements of optical dipole transition and study their behavior under axial 
electric field. We shall work with atomic units (a.u.), keeping for energy 1 a.u. = 
2Ry* = Ed [Ed is the effective ionization energy) and for distance 1 a.u. = a*^. 
A. Eigenenergies and eigenfunctions of HMI at zero bias field 
There exist several ways how to compute the eigenenergies and eigenfunctions of 
HMI. We consider one of them using variable separation in the time-independent 
Schrodinger equation followed by representation of solution in new variables via 
appropriate series expansions (see the work [29] and references therein). It is well 
known that the Schrodinger equation for HMI 

(V + E + - + -)vl/(r)=0, (37) 

where rA{B) denotes the distance between an electron and atom A{B), is sepa- 
rable in prolate spheroidal coordinates ^ = {rA + rB)/R, i] = {rA — rB)/R, ^ = 
arctan(y/x) (1 < ^ < oo, —1 <''7<1, 0<<y9<27r) and a (non-normalized) 
stationary electronic wave function of HMI can be written in the form \E' (^, 77, ip) = 
S (^) H (r^) $ {(f), where functions S (^) and H (77) meet the generalized radial and 
angular spheroidal wave equations, respectively, and the azimuthal function $ (ip) is 
the same as for hydrogen atom: $ {(p) = exp {im(p) j V^tt, where m = 0, ±1, ±2, .. is 
the azimuthal (or magnetic) quantum number. The most useful expansion for H (j]) 
is in a series of associated Legendre polynomials, whereas for S (^) the power series 
expansion is applied. Substituting these functions into corresponding wave equa- 
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tions and requiring them to be minimal solutions of these equations, we obtain the 
expansion coefficients and separation constants and that at fixed m depend 
continuously on the energy parameter E and on the internuclear distance R. To 
obtain electronic eigenenergies we now have to find the set of separation constants 
common for both radial and angular solutions for given m and R. Graphically, these 
can be sought as intersection points of curves from sets {C^{E)} and {C^(-E')} and 
corresponding values of parameter E will be the eigenenergies of HMI (for more 
details, see the Ref. [29]). 

The HMI eigenstates is completely specified by the set of quantum numbers 
{N^, Njj,m), where N^^N^j) equals to the number of zeros in radial (angular) func- 
tion S (H [r])). On the other hand, in the united atom limit i? — > the quantum 
numbers {nu,l,m) are used, where riu specifies the energy level, and / the angular 
momentum. Both sets are related by the formulas = 12^ — I — 1, = I — m. In 
what follows we shall use the second one for the state classification as more instruc- 
tive, i.e., {^,ri,ip) ^ \nulm). According to standard atomic notation, the letters 
s, p, d, /, g, h, ... are used to denote the values of / = 0, 1, 2, 3, ... in united atom 
limiting case, and the greek letters a, vr, 6, 0, ... denote the values of m. We can 
also classify the states of homonuclear HMI by the parity under the transformation 
r — r, namely, the state with symmetrical (antisymmetrical) wave function will 
be supplied with subscript g (u). 

Making use of computational framework sketched above, we have found the eigenen- 
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ergies and eigenfunctions of discrete part of HMI energy spectrum. At Fig. 2, the 
eigenenergies of first 20 a-states (m = 0) are plotted as functions of internuclear dis- 
tance i? for 20 < i? < 40 (on the whole, we found the eigenenergies of 64 states with 
171 = 0). We also computed the eigenenergies of states from subspaces with m 7^ 0, 
however, only a-states are relevant for the dynamical description of electron evolu- 
tion in the electric fields polarized along z axis that will be in the focus of further 
consideration. The HMI eigenenergies are found with absolute accuracy 10~^^ a.u. 
that is enough for this study, but the described algorithm enables one to calculate 
them with absolute accuracy 10^^*^ a.u. or even higher [29]. From Fig. 2 (a) it is 
seen that two low-lying states of HMI, i.e., \lsag) and |2pcr„), are almost degenerate, 
and the energy difference Ais^^^2pcTu = ^ (2pcr„) — e (Iscr^) decreases exponentially 
with R. The dependence of the electronic tunneling time tiso-<,-2p(t„ = ^/ 'Aiscrg-2pcru 
(here h = 6.582 x 10~^^ eVx s) between these states on R can be seen at Fig. 3, 
where the parameters E£)=-45 meV and a^=1.22 nm, taken from Ref. for Si:P2' 
in single valley approximation, are used. From this Figure we obtain an estimation 
on Tisag-2pau R ~30 um (:^20 a.u.) to be of order of 1 /is that is enough for 
carrying out some proof-of principle experiments, whereas for R ^50 nm (:^40 a.u.) 
one has for the tunneling time Tisag-2pau ~ 10^ s that would enable to run quantum 
algorithms. The states llscr^) and |2pcr„) are presented by even and odd superposi- 
tions of states localized on donors A and B. Inversely, the logical (localized) states 
can be expressed by even and odd superpositions of the states \lsag) and \2pau) of 
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HMI: |0(1)) = = i\lsag) ± |2pa„))/V2. 

The collection of low-lying electronic excited states falls into the subbands char- 
acterized at i? >>1 by the principal quantum number n of isolated hydrogen atom. 
The sets of HMI eigenstates, whose energies are pictured at Fig. 2 (b), namely, 
{|4/(T„) , \Mag) , |3pcr„) , \2sag)}, {|6/i(T„) , \5gag) , |5/cr„) , \4:dag) |4pcr„) , |3scTg)}, and 
{|7j(7u) , |6/i(Tg) , ..} (not labelled), correspond to n —2, n =3, and n —A, respec- 
tively. In their turn, these subbands can be subdivided further into the doublets 
\ksag) — |(A; -I- l)p(T„), /c=2,3,..; \kdag) — \{k + 1) fay), A;=3,4,..; and so on. The de- 
gree of hybridization/localization can be evaluated from state quantum numbers. 
For example, the pairs of states with I — s,p become degenerate at compara- 
tively small internuclear distances. The doublets of states with I — d, f are well 
resolved up to R ^30 and collapse quickly for larger R. We observe almost com- 
plete dissociation of HMI states, pertaining to the subband with n —2, into two 
pairs of states |2s^('^)) = {\2sag) ± |3pcT„) + \3dag) ± |4/ct„))/2 and |2p(T^(-^)^ 
(|2s(Tp) ± \Spau) — |3d(7p) ^ |4/(7u))/2 of isolated hydrogen atoms. Taking the aver- 
age over the energies of states pertaining to the subbbands with n = 1 and n = 2 and 
adding the energy of internuclear repulsion Snuci — ^/R, we arrive at isolated hydro- 
gen atom energies e{n) — 1/2'n?. On the other hand, the states with higher values 
of / remain non-degenerate on the whole interval of R under consideration. One 
can therefore treat those states as rehable transport channels to drive an electron 
between logical states. As we shall see below, the states \5gag) and \6hau) (whose 
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energies are plotted as thick red lines at Fig. 2 (b)) meet the conditions imposed 
on the choice of transport states by the high excitation selectivity and appropriate 
dipole coupling strength requirements. The states, whose energies are plotted as 
thick blue lines, may be exploited as transport ones either at R < 20 {\3dag) and 
|4/a„)) or at i? > 30 {\6gag), |7/ia„)), \7tag)). 

B. Matrix elements of optical dipole transition at zero bias field 
The numerical solution of Eq. (6) implies the knowledge of matrix elements {dkm} 
of optical dipole transition (ODT) between all states entering into the Eq. (4). 
In our simulation we restrict ourselves by consideration of electric fields polarized 
along z axis, so that the field-structure interaction Hamiltonian in Eq. (3) reads 
—eE {t) z. In this case, the matrix elements of ODT between states \k) and |m) can 
be calculated in spheroidal coordinates through the expression 



dkm = i^l^dv Jd^ [H, (r;) Ek (0] i-eR^vm [H^ (v) ^m (0] J (C, v) | , 

(38) 

1 oo 

Nkim) ^ Jdrjjd^ [lik(m) iv) ^Hm) (0]' J (e, V), J (e, ^) = - V', (39) 

-1 1 

where J {C,rj) is the Jacobian of the transformation from cartesian frame to spheroidal 
frame, and the relation z — (-R/2) ^rj is used. (The wave functions lik{m) iv) ^^'^ 
Sfc(m) of discrete part of HMI energy spectrum are real functions, hence the ODT 
matrix elements will be real as well). According to the selection rule Am = 0, im- 
posed by the structure axial symmetry and by the choice of the field polarization, 
and keeping in mind that logical states pertain to the cr-subspace of HMI eigenstates, 
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only the transitions among the states with m — are relevant. Besides, the analysis 
of Eq. (38) supplies us with other selection rule standing for allowed transitions Al 
to be an odd number. In other words, the states from subspaces {s, d, g, i, ..} with 
even I are dipole-coupled to the states from subspaces {p, /, h, j, ..} with odd I. It 
means that ODT selection rule upon Z, that follows from spatial symmetry of HMI, 
turns out to be relaxed in comparison with that of hydrogen-like atom for which 
Al = ±1. 

Using the equations (38) and (39), we have calculated the ODT matrix elements 
between all pairs of 64 low-lying cr-states of HMI as functions of internuclear distance 
R. In what follows, they will be used to define the right-hand side of Eq. (6). At 
the same time, in Eq. (21) we proceed with matrix elements between logical states 
|0) and |1) and excited ones. These matrix elements can be derived from expression 
do{i)k = [disagk^d2pcr^k) /V2. The dependencies of dok (1 a.u. = ea^) on R are 
given at Figs. 4 (a) - (c) for subbands with n =2, 3, 4. Note that in symmetric 
HMI dok — dik, if k denotes a state with even parity, and dofc — ~dik, if k denotes a 
state with odd parity. The calculation of matrix elements between localized states 
IgAB'^ from subband with n —1 (logical subspace) and localized states 
2p(7^'^) from subband with n =2, yields at i? = 38 dt^^2s ~ and df;\^ ^ 0.746 
that coincide with values dzis-2s = and dzis-2pa = {^o,*b) ~ 0.745 for hydrogen 
atom. 

C. Resonant and off-resonant electron dynamics in zero-bias case 



2s^'^) and 
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Here we present the results of numerical simulations on coherent electron evolu- 
tion, paying attention to performing ux qubit-state vector rotation. For simplicity, 
we work with square pulses and entering into Eq. (34) is the step function. To 
drive an electronic population in symmetric structure such as HMI, a single pulse is 
only needed, thus we integrate Eq. (6) with Hamiltonian H — Hq — eE^z cos {wot), 
and initial conditions corresponding to the localization of an electron into logical 
state |0), are Cuag (0) = C2pau (0) = l/ Our goal is to calculate the frequencies 
of Rabi oscillations between logical states and to estimate the degree of population 
leakage from logical subspace for different values of internuclear distance R, pulse 
strength Eq, and pulse frequency cuq- For this purpose, we analyze the probabilities 
Po (T) and pi (T) (po(i) {T) = Cis^g {T) ± C2po-„ {T) j 2) to find electron into the 
logical states |0) and together with total probabihty Ptr {T) = J2 Pk{T) of 
electron to be out of logical subspace, versus the pulse duration T. 

Let us introduce the dimensionless field energy e field = eEQa%/2Ry* and the 
dimensionless pulse duration T — 2Ry*t/h. According to the analysis given in Sec. 
II, the coupling coefficients that in atomic units take the form Ao(i)a; = Sfieiddo{i)k j 2, 
have to satisfy the requirements imposed by a concrete optical excitation regime. 
If one applies a resonant driving scheme, where HMI state |r) is used as transport 
state, the detuning of pulse frequency cuq from resonant frequency cuor must be 
much smaller than the value of coupling coefficient |Aor|- In its turn, |Aor| must be 
much smaller than LO^rw and order to minimize population leakage into the 
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states nearest-in-energy to transport state. Instead, the exploitation of strongly- 
detuned pulses (Raman scheme) implies the values of coupling coefficients |Aoa:| for 
all excited states to be much smaller than corresponding detunings. It is easy to see 
that one could attain this condition taking pulse strength as small as possible and 
detunings as large as possible. However, this brings about considerable reduction in 
Rabi frequency. In so far, we shall be interested in determination of optimal pulse 
parameter set, that would amount to rapid and robust implementation of quantum 
operations. 

We begin with resonant driving scheme that, being realized in symmetric struc- 
ture, may be considered as auxiliary one, since it is only able to inverse the popula- 
tion of logical states at discrete set of pulse durations when the electron is concen- 
trated into logical subspace. The electronic resonant population transfer in three- 
level and in four- level structures was studied theoretically in Refs. [21] and [23] , 
where the probabilities Po, Pi, and ptr were found at exact resonance (zero detuning 
from transport level) to be 

PO (T) = cos^ , pi (T) = sin^ , p. (T) = 1 sin^ . 

(40) 

One can observe from Fig. 5 that these expressions are in excellent agreement with 
numerical curves. We have performed our simulations on resonant dynamics for sev- 
eral values of e field (different pulse strengths Eq) and pulse frequencies uq, matching 
resonant frequencies of HMI (different dotr) in order to define the dependency of Rabi 
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frequency on these parameters. We have revealed that the pulses with Sfieia < 0.003 
being tuned on resonance with states \5gag) or \6hay), provide good selectivity, low 
population leakage at the end of NOT operation (pi {Tj^qt — '^/'^^r) > 0.999 ), 
and quite high speed of electronic transfer {T^qt ~ 10^ that is of order of hundreds 
of picoseconds for Si:P2" structure) .We have also performed numerical simulations 
on the electron dynamics for R — 22 in order to check the possibihty of using 
the states from HMI subband with n = 2 as transport ones. As it was expected, 
the resonant population transfer between logical states via excitation of one of the 
states {\'idag) , \Afau)} do really take place. However, large dipole moments for 
those transitions (see Fig. 4 (a)) require the field energy to be sufficiently low since 
the energy spacings between subband levels are still small. As a result, the time 
needed for complete electron transfer is of the same order as it was for larger R, 
where the states from third and fourth subbands, characterized by smaller values of 
ODT matrix elements but demonstrating higher resolution over energy, play role of 
transport channels. Thus the formula VIr — |£'odotr|/4 = |Aotr|/2 approximates with 
high accuracy the Rabi frequency for resonant electronic transfer for 20 < i? < 40. 

Now we present the results of calculations for the off-resonant scheme. Typical 
curves for po, Pi, and ptr, reflecting coherent electron dynamics driven by strongly 
detuned pulses, are plotted at Fig. 6. They demonstrate essentially two-level oscil- 
latory behavior where the population is mainly localized in logical subspace. This 
picture differs from that obtained for resonant scenario by sharp decrease in the 
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amplitude of population combined with considerable increase in the oscillation 
period of po and pi (by order of magnitude or more). According to Eqs. (21), (31), 
(33) and (34), in the symmetric structure Aq — Ai = 0, ©o — 7''/2, and 

po (T) ^ cos' (QrT) , pi (T) ^ sin' (fi^T) , (41) 

where — IA2I (see Eq. (25)). Relative to the transport states, we are only able 
to estimate under adiabatic elimination an order of magnitude of total probability 

as ptr ~ (Aofc/^ofc)Lax << ^ its oscillation frequency fltr ~ (5ofc/Aofc)max » 
so that Ptr exhibits fast oscillations with small amplitude. 

In order to compare the Rabi frequencies found numerically with those calculated 
within analytical framework, we plot both data types on Fig. 7 versus pulse fre- 
quency (shifted by £0) for R — 30 and R ~ 38 and for the field energies Sfieu— 0.003 
and e field— 0.005. To extract numerical values of IA2I, we fit the curves for pi (T) by 
squared sine function and then define as the sine frequency. These values, drawn 
as filled circles, correlate well with analytical results of Eq. (21) pictured by solid 
curves. Since the approximation applied in derivation of Eqs. (20) and (21) do not 
allow one to work in near neighborhood of HMI energy levels, we have left empty the 
intervals around the levels marking them by vertical dotted fines. We shall regard 
the oscillations as two- level ones if the depth of modulations of po and pi, arising 
from the non-resonant population of excited states, is smaller than 0.01. At Fig. 7, 
there are several points (enclosed by open red squares) satisfying to this conventional 
criterion. Other numerical data points correspond to the ptr maxima (or modula- 
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tion depths) ranging from 0.01 to 0.05. For the points located at the boundaries 
of the intervals, on which approximate solution is valid, the maximum values of ptr 
turn out to be 0.05 or higher. With further approaching of the pulse frequency to 
one of HMI resonant frequencies, the off-resonant oscillatory picture transforms into 
resonant one. The off-resonant dynamics aX R — 22 (not shown) is quite similar 
to that presented at Fig. 7. Again, almost ideal two-level oscillations are obtained 
if we tune the pulse frequency into the middle of doublet {{bgag) , \Qhau)}- With 
that, the use of doublet {\?>dag) , |4/(T„)} becomes possible provided that the pulse 
strength is not high {s field < 0.001). As it will be shown below, we cannot efficiently 
operate with such pulses since the Rabi frequency of qubit rotations appears to be 
very small. 

Note that the character of oscillations of total population of excited states is regu- 
lar enough in both resonant and off-resonant cases; this observation can be explained 
by assumption that only several excited states participate the dynamics. Such an ex- 
planation is obvious for resonant excitation scheme, but in off-resonant case, where 
all excited states are equivalent in dynamical sense, rigorous arguments are needed. 
To find them, let us consider the sum A2 in more details. As we have mentioned 
before, the excited states of HMI can be classified by tunnel coupling strength be- 
tween symmetric and antisymmetric states pertaining to the same doublet. For 
weakly hybridized states \ki) and \k2), provided that the value of tunnel splitting is 
much smaller than pulse detuning from one of these states, we can write ~ 6^2- 
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We suppose, that the pulse energy bandwidth is much larger than \5k^ — 5k^\ so we 
are disable to resolve the doublet states under consideration. Furthermore, in this 
case Xoki ~ Aofca, '^ofci ~ ^ikn and Aofcj ~ —^ik2 (see Sec. IV B). Therefore, the 
terms Xok^ik/^k arising from these states have opposite signs and cancel each other, 
so that the contribution of the doublet , \ k2)} to A2 is minimal. To estimate 
the contributions from higher states with energies e > —0.05, we should take into 
account two circumstances, i) The ODT matrix elements between logical states and 
excited states with equal I decrease with the energy growth whereas the detunings 
increase with the energy growth, ii) High-lying states are closely spaced to each 
other and the distance between two neighboring states (not necessarily from same 
doublet) decrease rapidly with the energy growth. Thus, one may expect that cor- 
responding terms in A2 will either cancel out each other, as it was for degenerate 
states, or have insignificant effect on the sum convergence because of their subse- 
quent reduction versus state energy. We think that both issues are important and 
their cooperative effect takes place, establishing our observation on the oscillation 
type of ptr- Strictly speaking, for long times T, when flnT » 1, this is not true, 
because an internal structure of the sum ptr, containing the oscillating terms of 
different but close frequencies, causes these oscillations to be averaged out and ptr 
tends to its average value ~ max (pt^) /2. 

Utilizing similar arguments, we specify an algorithm of choice of the pulse fre- 
quency and amplitude in off-resonant driving regime. First of all, we determine 
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the doublet whose states kij and /C2) are characterized by sufficiently large ODT 
matrix elements and high tunneling rates (large value of the energy splitting in 
comparison with the pulse energy bandwidth). Then, by adjusting the pulse fre- 
quency to the center of energy gap between doublet states, we obtain for detunings 
S^^ — — ^^^^ case, the terms Xok^ik/^k enter into A2 with the same sign 
resulting in amplification in the sum. Doing so, we should remember that the field 
energy s field (i-e., pulse strength) must be much smaller than the doublet energy 
splitting. As it can be seen from Fig. 2 (b), upper bound value for Efieia is to be 
set ~ 0.005 at the whole considered range of R. Actually, it means that we can reli- 
ably operate with frequencies that fall within the interval containing the states from 
HMI subbands with n — 3 and n — 4. Moreover, the frequency choice just outlined 
cannot guarantee low population of excited states. A general way to minimize this 
population implies the use of moderate pulse strengths at fixed detunings. Note, 
that two-level Rabi oscillation regime with max(pjr) < 0.01 can be achieved under 
this frequency choice, since marked points on Fig. 7 lie in the middles of energy 
gapes separating the doublet states \5gag) - |6/i(7„), \6gcrg) - \7hau), and \7iag) - 
|8j(7„). With that, the pulses have to be strong enough in order to perform qubit 
rotations in the times that are shorter than decoherence time. For the off-resonant 
driving scheme, where the decoherence effects due to the finite population of excited 
states (e.g., the electron- phonon relaxation) are reduced considerably, the dephas- 
ing of localized states is expected to be the main source of coherence losses. In our 
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simulations we require that the operation times (Top ~ do not exceed 10~^ s. 

It gives a lower bound value on the pulse strength to be ejieid ~ 0.001. For SiiPa" 
structure parameters, the allowed range 0.001 < e field < 0.005 corresponds to the 
pulse strengths of hundreds of V/cm that is by two orders smaller than the field 
strengths considered in the voltage-based driving schemes [12] - [16]. 

Provided that the pulse parameters are selected correctly, the probability ampli- 
tudes Co and ci (and, consequently, probabilities po and pi) evolve smoothly against 
time thus allowing one to continuously rotate qubit state vector along fixed meridian 
with if = it/ 2 on the Bloch sphere, namely, to create the superposition of logical 
states of the form |^ (t)) = ^po (t) |0) -i^Jpi (t) |1), where po (t) + pi (t) ^ 1. Con- 
trol over relative phase between logical states, needed for creation of an arbitrary 
superposition state of the qubit, requires the logical states to be addressed indepen- 
dently by two different pulses to maintain Aq 7^ Ai in Eq. (20). To achieve this goal, 
one should break the central symmetry of HMI making use of external field, thus 
the important quantum operations given by Eqs. (34) - (36) will be realized in this 
case after application of bichromatic pulse. Besides, there exists a more familiar way 
utilizing the adiabatically varied bias field to produce a shift between the energies 
of logical states. 

D. Phase and population dynamics in nonzero-bias case 
As it was demonstrated in the works [5], [IS], [2D], and [3D], an uniform elec- 
trostatic field, polarized along z axis, disturbs spatial symmetry of DD structure 
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relative to the origin and produces a necessary energy shifts. The Hamiltonian of 
HMI, subjected to the action of an external axial field, reads H = Hq — cEqz. We 
compute the eigenfunctions | ^i)} of HMI with the field, expanding them over the 



eigenfunctions without the field: ^i) — Z^Cfei |^fc)- The eigenenergies and 

' k 

the expansion coefficients Cki are calculated directly through the diagonalization of 
the Hamiltonian in the basis of 64 u-states, found in Sec. IV A. Similarly, the ODT 
matrix elements between the HMI eigenstates in the non-zero bias case are 

expressed via the ODT matrix elements between the HMI eigenstates in the zero 
bias case, {dmn}, as follows: dij = E C^iCnjdmn- 

m,n 

The dependencies of eigenenergies on static field energy ef^^^i^ (defined in com- 
plete analogy with e field) are given at Fig. 8 for internuclear distance R — 38. 
With respect to the hybridization degree of HMI states at efi^id — 0, the important 
difference in energy behavior of states, originated from localized and delocalized 
states, is observed at £^*eZd 0- that the energies of states, corresponding 

to weakly hybridized states of HMI without bias field (e.g., the states from sub- 
bands with n = 1 and n — 2), are approximately linear functions of sfi^id {due 
to Stark effect). The generation of the energy difference A — e-i — Eq = efi^i^R, 
where R — 38, requires an application of bias field with the strength Eq 10 V/cm 
(see Fig. 8 (a)). If ef^^id > 0' states, whose energies bring down (up) with 
the static field strength, are locahzed predominantly on "deep" donor B ("shallow" 
donor A). From Fig. 8 (b) one can observe a general tendency of the HMI spec- 
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trum to become complicated with the field growth due to both the full lifting of 
Coulomb degeneracy (in contrast to the partial lifting in HMI spectrum without 
field) and the hybridization of states pertaining to different subbands (giving rise 
to crossing-anticrossing pattern on Fig. 8 (b)). Otherwise, the energies of delocal- 
ized states |5(7cr) and \Qha) (thick red curves on Fig. 8(b)) vary with static field 



energy as e [bg {6h) a] ^ l[e {5gag) + e (6/i(t„)] t |y {^fieidR) + ^Iga.-sha^ , where 
^5gag-6hau = ^ {QhcTu) — s {^gcTg) is the tunucl splitting at zero bias field. It means 



Afield 



< 0.001. 



that the tunneling rate for those states remains sufficiently high if 
Therefore, we shall continue to operate with those states exploiting them as the 
transport channels. The dependencies of the matrix elements on the field strength 
for transitions, connecting logical states and transport states l^gcx) and \6ha), are 
given on Fig. 9. The examination of the values of matrix elements confirms the fact 
that the maximum of electron density for the state l^ga), originated from symmetric 
state \5gag), is displaced onto "deep" donor B, whereas for the state \6ha), origi- 
nated from asymmetric state \6hau), it goes onto "shallow" donor A. The deeps on 
the curves appear due to the electron density redistribution at crossing/ anticrossing 
points, where transport states interact with localized states (see Fig. 8 (b)). 

In Ref. [T7], it was demonstrated that elementary single-qubit operations can be 
performed in asymmetric DD structure, driven by two resonant pulses. Here we 
check this proposal for zero-detuning case by the simulation of NOT operation. Ac- 
cording to our previous results, the resonant asymmetric three-level scheme requires 
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that \X(ytr\ — \Xitr\, \Eo/ Ei\ — \ditr/dotr\- Since the ODT matrix elements for dis- 
tinct arms of the excitation scheme differ from each other, we need to compensate 
the difference in couphng coefficients Xotr and Ai^^. by choosing the pulse strengths 
so that to fulfill above condition. Now we integrate Eq. (6) with the right-hand side 
defined by the set | of the ODT matrix elements and with the initial conditions 
reflecting the localization of electron at the beginning of pulse action in the ground 
state |0). As usual, the probabilities po{T), pi{T), and ptr{T) (Fig. 10) illustrate 
the qubit state inversion for the case when both pulses are in the exact resonance 
with the state \5ga) (the other parameters are indicated at the plot). We can re- 
veal that the arms of excitation scheme remain inequivalent, despite of that the 
couplings have been balanced. The time dependency of pi is smooth, but the plot 
for pq demonstrates a fringes superimposed on ideal "three-level" curve (compare 
with that pictured on Fig. 5) and originated from non- resonant excitation of "deep" 
donor states whose energies lie closely to the \5ga) states. For £^*g;^=0.0004, the 
difference between energy of transport state and those of the states of donor B is 
~ 0.006, whereas the energies of states localized on donor A lie by ~ 0.02 higher 
than the transport state energy. Add, that the character of dynamical picture is not 
changed noticeably, if one uses the state \6ha). To suppress these fringes and refine 
oscillation picture, one should apply the pulses with lower strengths. In spite of 
presence of these perturbations, the electronic transfer probability is very high and 
the period of the oscillations is short enough. We expect this type of evolution to be 
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conserved for various quantum operations carried out in asymmetric DD structure 
under the resonant excitation regime. 

Unlike in the case of symmetric HMI, here we were unable to find pulse parameters 
for asymmetric off-resonant driving scheme that would provide a robust and fast 
implementation of basic quantum operations. This lack takes place because the 
pulse with frequency cuo, being tuned slightly below the state \5ga) for the transition 
|0) \5ga), addresses unwanted nearly-resonant transitions connecting the state 
|1) with the states localized on "shallow" donor. Thus, the selectivity requirement 
is violated. The attempts to achieve the reliable two-level oscillation picture have 
forced us both to reduce the pulse strengths and to vary the detuning in such a 
way that both pulses would drive only the transitions prescribed by ideal excitation 
scheme. The goal is likely to be attained by this method, but the characteristic times 
will be inappropriately long in comparison with decoherence time. In particular, the 
asymmetric off-resonant scheme, described in Sec. Ill, will function correctly if A is 
much smaller than the energy spacings between the states lying around the energy 
So + ujq. At same time, the coupling coefficients must be much smaller than A, 
in order to address separately each arm of driving scheme by its own pulse. Prom 
the Fig. 8 (b) we obtain the estimations ef^^^i^ < 10~^ and A < 0.004 that, in its 
turn, implies |Ao(i)fc| to be less than 10^^. It tells us that under these conditions 
the off-resonant Rabi frequency Qr is of order of 10^ (or smaller) that is out of 
interest of this paper. 
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Actually, it means that we can only deal with symmetric DD structure where the 
possibility of realization of the off-resonant qubit-state ax rotation has been already 
demonstrated. However, we know this type of quantum evolution is not sufficient 
for implementation of arbitrary quantum operation, since the phase control is also 
required. This difficulty may be overcome with the help of direct introduction of the 
phase difference between logical states by the voltage pulse that produces necessary 
energy shift. This technique is widely used in quantum dot structures to align the 
energy levels of different dots that is needed for experimental investigations of their 
transport and spectral properties. Here we consider a simplified model of voltage 
pulse action on the HMI replacing non-isotropic electric field, generated by the gate 
Va (Fig. 1 (a)), by uniform axial electric field, as it is demonstrated at Fig. 1 
(b). The time dependency of voltage pulse is approximated by the step function. 
In this case, elementary calculations predict the qubit phase oscillations with the 
frequency O^p^ase = A/h. The general expression for qubit-state evolution operator 
in the two-level approximation, accounting for different voltage pulse shapes, can 
be found in Ref. [30|. We illustrate this type of quantum dynamics at Fig. 11, 
where two-level Rabi oscillations between the states 1+)^ = (|0) + ^ y/2 and 
\^)h ~ ~ obtained from logical states by the Hadamard rotation, are 

presented. Relative to computational basis, these oscillations indicate on the phase 
dynamics generated by the operator az- For example, the application of voltage 
pulse during one-half oscillation period amounts to the phase shift by vr between the 
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states |0) and |1) and takes about several picoseconds. With that, the population 
leakage from logical subspace remains negligible {ptr < 10^^). At the intersection 
points, the Hadamard superpositions collapse into one of logical states. 

Summarizing the results obtained in this and previous subsections, we conclude 
the numerical study of electron dynamics by observation that the strategy, utilizing 
sequential implementation of ax and az rotations, appears to be more reliable than 
the strategy based on simultaneous action of voltage and optical pulses. With that, 
one may regard the formulas deduced in Sec. Ill for asymmetric off-resonant scheme 
as describing non-ideal nearly-symmetric scheme driven by single pulse, where other 
external fields (e.g., acting on nearest qubits) slightly disturb the spatial symmetry 
of structure in uncontrollable way. In this case, small differences in frequencies of 
the arms quantified by the asymmetry parameter A, will amount to quantum error. 

Of course, we should take into account that the driving scheme of Eqs. (34) - 
(36), whose efficiency has not been confirmed above for effective HMI, can describe 
the qubit dynamics for other models of DD structure (for example, that of heteropo- 
lar molecular ions [16]). The main issue, that makes the verification of proposed 
algorithm of quantum control difficult, is concerned with necessity to know in de- 
tails the energy spectrum and the ODT matrix elements for each structure under 
consideration. 

V. DISCUSSION 

To provide more clarity in the understanding of the advantages of the charge qubit- 
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state engineering presented above, let us compare the resonant and off-resonant ex- 
citation scliemes. Tlie resonant optical driving of the DD structure modelled by 
single three-level scheme has been studied in Refs. [T71 HB] . In the case considered 
here a more complex dynamics can take place involving more than one three-level 
scheme. For example, if we tune the lasers on resonance with the transition between 
the logical subspace {|0) , |1)} and a state |r) , r G {k}^^^ located near the top of 
the barrier, a number of states with the energies close to Er will be excited as well. 
This picture is quite expected in the hydrogen-like molecular ions for the high-lying 
states which energies are within the interval Asr < Ao(i)r • However, for the sym- 
metric structure it means that there will be no transitions between the qubit states. 
It is because the exited states belonging to the same doublet are presented by the 
symmetric and antisymmetric superpositions of the excited states of isolated donors 
which, being excited simultaneously, interfere constructively on one donor and de- 
structively on another one. As a consequence, the donors are excited independently. 
This effect becomes more significant as the interdonor distance R increases and the 
tunnel coupling between the donors decreases. We have observed such type of elec- 
tronic dynamics when pulse was tuned into resonance with low-lying excited states 
with / = s,p,d,f. When the energy splitting of the maximally-resolved doublet 
becomes comparable with the coupling coefficients of the optical dipole transitions, 
the process of the electron transfer between the donors is terminated. Note that 
the optically driven DD structure will demonstrate the similar behavior if one of 
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the pulses is short enough so that its duration is T < l/Aerr and thus it contains 
harmonics in the frequency range ^cj ~ 1/T > Ae^. Again, the states with the 
energies belonging to the interval Ae^ will be excited simultaneously giving rise to 
the electron transfer blockage just outlined. In our simulations, we have arrived at 
this regime via application to HMI a short and intense {e field ~ 0.1) pulse. 

The reliable resonant scheme thus deals with single transport state (for HMI, the 
states 5(7cr(g)^ and 6/icr(u)^) and is very sensitive to the pulse detuning from the 
resonance with that state. For example, the non-zero detuning always produces 
an amplitude error in NOT gate because of incomplete depopulation of the initial 
state when the pulse is off [21j. On the contrary, the use of the off-resonant pulses 
enables one to exploit the whole number of excited states (from which only several 
ones participate substantially in electronic dynamics). Moreover, we don't need 
to control the pulse frequencies with high accuracy since a small variation in the 
pulse detunings brings about an insignificant change in the Rabi frequency (see Fig. 
7). The computational errors originated from the frequency renormalization can be 
corrected by the corresponding change in the pulse duration due to the smooth time 
dependencies of the probability amplitudes, as it was shown at Fig. 6. The only 
requirement that must be followed closely for successful electron state manipulations 
is the Raman two-photon resonant condition (9). 

The selectivity of the electron resonant transfer requires also a strict control over 
the pulse polarizations. The transport states in the molecular ion are formed through 
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the hybridization of those individual donor states whose wave functions are extended 
along the axis z that coincides with the interdonor direction. Other states (e.g., vr- 
states) are hybridized weakly and cannot assist efficiently in the electron dynamics. 
Their excitations are due to the pulse components polarized along the axes x and 
y. It amounts to the population leakage into the non-hybridized single-donor states 
with the energies lying in the close proximity to the energy of the transport state. 
Let us define the small angles 7„x and 7„j, that characterize the deviations of the 
ra-th pulse polarization from the axis z: 

E„, = E„2 + E„2; cos (7r/2 + 7„^) + Y^^y cos (7r/2 + 7„y) , , . 

\lnx\ , \lny\ < 1, n = 0, 1, 

then the probability of successful implementation of the quantum operations is re- 
duced by a factor of w ~ 1 — max(7^^, 7^^). In the off-resonant case, the populations 
of those states remain negligibly small (~ \\nk\ j ^X) ^-nd the corresponding channel 
of population leakage is blocked. 

The important difference between the resonant and off-resonant schemes lies in 
the treatment of the decoherence problem. We know the relaxation rates from 
the transport state caused by the spontaneous photon/phonon emission during the 
resonant excitation [TS] may be high enough to corrupt the qubit state. In the 
off-resonant scheme the population of the intermediate state(s) is negligible and 
the probability of relaxation is drastically reduced. The influence of the residual 
population of the intermediate state on the adiabatic electron transfer in the three- 
level scheme was examined in Ref. [31] for the gaussian pulses. It was shown that 
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the error introduced by the spontaneous emission together with the error due to the 
non-adiabaticity are inversely proportional to the pulse detuning and can be made 
small enough to allow the fault tolerant quantum computation. 

Note that the complete population transfer between the qubit states, or NOT op- 
eration, requires that Aq = Ai. This is naturally met for nearly symmetric DD struc- 
tures where A ^ and |doA:| ~ |dife|. In general, however, one should keep in mind 
that |dofc| 7^ \dik\ that makes the performing of the condition Aq = Ai very prob- 
lematic. It seems then reasonable to point the other way for the population transfer 
based upon the pulse-shaped techniques. Such methods, e.g., the stimulated Raman 
adiabatic passage (STIRAP) [26], are very robust against the pulse/structure im- 
perfections and would allow one to handle the quantum information carefully. The 
theory of the adiabatic population transfer via multiple intermediate states, includ- 
ing the off-resonant case, was presented in Ref. [27] . Note that for the pulses strongly 
detuned from the resonance, the time ordering is no more important since successful 
population transfer may be attained for both intuitive and counterintuitive pulse 
sequences. If initially cq (to) = 1, ci (to) = 0, the intuitive (counterintuitive) pulse 
ordering means that lim [fo (t)/ fi (t)] = oo (0) and lim [/o (t)//i (t)] = (oo) and, as 
it follows from Eq. (24), (to) = (vr), 6 (T) = vr (0). The population transfer may 
be understood as the adiabatic temporal development of the eigenstate |-|-) (|— )) 
for the intuitive (counterintuitive) pulse ordering. As it is seen from Eq. (32), the 
qubit state inversion is realized in the asymmetric DD structures if the conditions 
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arg [A2 (T)] ± Cl (T) = 7m and TA = tt (2m + 1) are fulfilled. The detailed analysis 
concerning the arrangement of the pulse shapes in STIRAP can be found elsewhere 

The effect of the detuning 5two-ph from two-photon resonance (the spacing between 
dashed horizontal lines on Fig. 1 (b)) should be also taken into account at more 
profound level of investigations. In nearly-symmetric off-resonant case (|Ao — Ai| <^ 
IA2I) the problem enables analytical solution. If electron was initially localized into 
the state |0), the total excitation probability (including transfer probability into state 
|1)) at large detunings (IA2I <^ \5two-ph\) is of order of {1^2/ ^two-phf ^ 1- Therefore, 
if the detunings from two-photon resonance are significant, electronic excitations 
from localized state |0) do not occur. In the opposite case of small detunings, i.e., 
when \5two-ph\ ^ 1^2 1, the probabihty of successful electronic transfer is given by 
formula ^ 1 — {5two-ph/'^^2f' (provided that other parameters are chosen in 
optimal way). 

Finally, in our analytical treatment we consider continuum states only as virtually 
excited transport channels, neglecting the possibility of two-photon resonant elec- 
tron transitions from logical subspace into the continuum states, that reside within 
interval around the energy £21^0 = ^0 + 20^0, via intermediate low-lying excited states 
(not necessarily highly populated) with energies ~ £0 + ^^o- It is known that in 
isolated hydrogen atom the matrix elements for the transition ^ |2j9cr) and the 
matrix elements for transitions connecting the state \2pa) with s (or d) continuum 
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waves with energies pertaining to the interval pointed above, have the same order of 
magnitude. Hence, one could expect the electron dynamics in DD structure induced 
by optical driving pulses to be much more complex. However, in experimental in- 
vestigations of low-frequency excitations in atomic systems these processes do not 
manifest themselves at the level that would establish the necessity to include the 
continuum into theoretical model. Perhaps, it could be explained by the arguments 
we have used in Sec. IV C to account for our numerical results reflecting very small 
influence of high-lying states of discrete part of HMI spectrum on electron dynam- 
ics. In principle, the continuum states being addressed directly can also be used as 
transport channels, but, at the same time, can bring about additional decoherence 
(see, e.g., [32] )■ 

VI. CONCLUSIONS 

In this paper we have considered the one-electron double-donor structure sub- 
jected to the action of optical and electrostatic pulses. Unlike the other systems 
proposed to serve as the potential candidates for the solid state optically-controlled 
qubits (double quantum dots, rf-SQUIDs), the double-donor structure is character- 
ized by sufficiently high density of the bound states at the edge of the barrier that 
separates the donors. It means that the three-level resonant scheme proposed ear- 
lier to implement the desired qubit-state evolution may be unsuitable to maintain 
the appropriate selectivity of the optical excitations. On the other hand, the off- 
resonant scheme looks as more efficient for the qubit manipulations and robust in 
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comparison with the resonant scheme. Though the Raman evolution of the qubit 
is slower than that in the case of the resonant driving, it seems to be more reli- 
able for the implementation of quantum operations. We have shown that the basic 
single-qubit operations may be performed on the DD structure for several pulse 
and structure parameter choices. Numerical simulations, carried out on structure 
modelled by effective hydrogen molecular ion, have confirmed the validity of our an- 
alytical framework where three important simplifications, concerning the treatment 
of dynamical problem, have been made. There are: i) rotating-wave approxima- 
tion, ii) adiabatic elimination, and iii) neglecting the transitions between high-lying 
states. Although all of them are widely used, it is, to our knowledge, for the first 
time when these approximations are verified in rather complicated dynamical study 
of multilevel system. 

The information about the structure and pulse parameters is contained in the 
Rabi frequency of the two-level oscillations. This frequency can be defined experi- 
mentally for each set of the detunings, the strengths, and the durations of the pulses. 
The results of those measurements could be used to reconstruct the features of the 
spectrum of the DD structure. 

Note that the method of the electron-state manipulations by optical means can 
be applied also to the spin-based encoding schemes like that of Ref. The 
implementation of optically controlled effective electron spin exchange described in 
Ref. |E3] for the two-electron double-dot structure, can be generalized on the two- 
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electron DD structure. The use of the adiabatic schemes is of the particular interest. 
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FIGURE CAPTIONS 

Fig. 1 (color online), a) Schematics of the quantum state manipulation in the 
one-electron double-donor structure. A pair of (phosphorous) donors A and B (one 
of them being singly-ionized) implanted into semiconductor (Si) matrix is addressed 
by optical pulse (s). Additionally, electrostatic pulse generated by voltage gate Va 
varies the structural potential. The desired final orbital state of electron is attained 
due to cooperative effect of both pulses on the structure, b) Potential profile of 
DD structure, modelled by effective hydrogen molecular ion (see Sec. IV), along 
structural axis z. The qubit states |0) and |1) are defined by the localized orbital 
states of the donors A and B with the energies Eq and £i, respectively. They are 
coupled to a collection of excited states by two optical pulses with the frequencies 
cuq and cui. Uniform axial electrostatic field is applied across the structure in order 
to break the symmetry and to introduce the energy difference A = £i — £o- Here, 
the distance R between donor centers is 38 a.u. and the electrostatic field energy 
afield = 4 X 10"^ a.u. so that A = x £^*g,^ = 0.0152 a.u. The energies of the whole 
molecular ion are obtained from electronic eigenenergies -|- energy of internuclear 
repulsion 1/R. 

Fig. 2 (color online). The dependencies of electronic eigenenergies of 20 low-lying 
eigenstates from cx-subspace (m = 0) of hydrogen molecular ion on the internuclear 
distance R. a) Ground-state energies e{lsag) and e{2pau). b) Excited-state ener- 
gies. Here, the eigenenergies of transport states \5gag) and |6/i(T„) are presented by 
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thick red curves whereas the eigenenergies of the states which may also be used as 
transport ones are drawn by thick blue lines. The vertical dotted lines correspond 
to the internuclear separations R— 22 a.u., R — 30 a.u., and it! = 38 a.u., for which 
dynamical simulations are performed. 

Fig. 3 (color online). Characteristic electronic tunneling time Tisag-2pay, between 
the localized states \lsag) and \2pau) of hydrogen molecular ion, calculated for 
parameters that correspond to single- valley approximation for Si:P^ double-donor 
structure, as function of the internuclear distance R. 

Fig. 4 (color online). The values of matrix elements of optical dipole transitions, 
connecting the logical state |0) and the low-lying excited states pertaining to the 
hydrogen molecular ion subbands with a) n = 2, b) n = 3, and c) n = 4, as functions 
of the internuclear distance R. 

Fig. 5 (color online). An example illustrating the resonant qubit-state manipula- 
tion. The time dependencies of probabilities Po{T) and Pi{T) to find electron into 
logical states |0) and |1) and that of total probability Ptr{T) to find electron into 
excited states are plotted for initial condition co(0) = 1, Cjk^o(O) = 0. The numerical 
solutions are given by solid curves, and the analytical solutions, presented by equa- 
tion (40) and supplied with primes at the insets, are visualized by dashed curves. 
The pulse is in exact resonance with transport state \6hau)- Other parameters are 
given at the plot. 

Fig. 6 (color online). An example illustrating the off-resonant qubit-state manipu- 



50 



lation. The probabilities Po{T), Pi{T), and Ptr{T) demonstrate almost ideal two-level 
oscillation picture. The maximum of total probability of population leakage from 
computational subspace is 0.01. The pulse frequency is tuned into the middle of 
HMI subband with n — 4 (red-squared point on Fig. 7 (b)). Other parameters are 
given at the plot. 

Fig. 7 (color onhne). The Rabi frequency IA2I for the off-resonant symmetric driv- 
ing scheme is plotted vs the pulse frequency lj = ujq at two values of pulse strength 
Afield — 0.003 a.u (blue), and S/iejd = 0.005 a.u. (green) for internuclear distances 
a) = 30 a.u. and h) R = 38 a.u. The curves correspond to analytical results of 
Eq. (21) whereas full circles mark the values of IA2I extracted from numerical data. 
Rabi frequencies, at which the probability ptr of leakage from logical subspace is 
lower than 0.01, are enclosed into red open squares. Vertical dotted hues designate 
the energy levels near which off-resonant approximation becomes inapplicable. 

Fig. 8 (color online). The electronic eigenenergies of HMI plotted in the depen- 
dence on static field energy £j*eZd -R = 38 a.u. Dotted vertical line designates 
the value sf^^^i^ = 4 x 10~^ a.u. for which the dynamical simulation is carried out. 
a) The energies £0 and £1 of logical states, b) The energies of excited states. Thick 
red curves denote the eigenenergies of doublet states \5ga) and \6ha) that remain 
delocalized in the presence of electrostatic field, and can be exploited as transport 
channels for resonant manipulations on electronic wave function. 

Fig. 9 (color online). The matrix elements of optical dipole transitions, connecting 
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the logical states |0) (dashed curves) and |1) (sohd curves) with the states \5ga) and 
\6ha), in the dependence on static field energy sf^^^i^ > for = 38 a.u. 

Fig. 10 (color online). Resonant manipulation of electron orbital state in asym- 
metric HMI at the electrostatic field energy £j*e/<i = 4 x 10""^ a.u. Both pulses are 
in exact two-photon resonance with state \5ga) and their amplitudes (field energies 
£ field Q and e field i) are chosen so as to equalize the coupling coefficients of effective 
three-level excitation scheme (see text for details). 

Fig. 11 (color onhne). The populations p+{T) and P-{T) of Hadamard-rotated 
states 1+)^ and vs electrostatic pulse duration T. Two-state Rabi oscillations 
in the Hadamard frame illustrate dynamics of relative phase between logical states 
|0) and |1). The population leakage into excited states is negligible. 
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Fig. 1 (b) 
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